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Abstract: Water resources are precious in arid and semi-arid areas such as the Wadis of Iran. To sustainably 
manage these limited water resources, the residents of the Iranian Wadis have been traditionally using several 
water use systems (WUSs) which affect natural hydrological processes. In this study, WUSs and soil and 
water conservation measures (SWCMs) were integrated in a hydrological model of the Halilrood Basin in 
Iran. The Soil and Water Assessment Tool (SWAT) model was used to simulate the hydrological processes 
between 1993 and 2009 at daily time scale. To assess the importance of WUSs and SWCMs, we compared 
a model setup without WUSs and SWCMs (Default model) with a model setup with WUSs and SWCMs 
(WUS-SWCM model). When compared to the observed daily stream flow, the number of acceptable 
calibration runs as defined by the performance thresholds (Nash-Sutcliffe efficiency (NSE)20.68, —25%<percent 
bias (PBIAS)S$25% and ratio of standard deviation (RSR)<0.56) is 177 for the Default model and 1945 for 
the WUS-SWCM model. Also, the average Kling—Gupta efficiency (KGE) of acceptable calibration runs 
for the WUS-SWCM model is higher in both calibration and validation periods. When WUSs and SWCMs 
are implemented, surface runoff (between 30% and 99%) and water yield (between 0 and 18%) decreased 
in all sub-basins. Moreover, SWCMs lead to a higher contribution of groundwater flow to the channel and 
compensate for the extracted water by WUSs from the shallow aquifer. In summary, implementing WUSs 
and SWCMs in the SWAT model enhances model plausibility significantly. 
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1 Introduction 


Over-exploitation of water resources and water scarcity have become a worldwide and prevalent 
problem in most arid and semi-arid regions such as the Wadis in Iran (Voss et al., 2013). In addition 
to a strong variability of precipitation in space and time in arid and hyper arid climates, these 
regions experience the lowest amounts of precipitation in Iran (Khalili and Bazrafshan, 2004; 
Motiee et al., 2006). In arid environments, most of the rainwater is lost to evaporation (Tavakoli et 
al., 2010). Under such circumstances, the limited available freshwater is a fundamental and 
valuable resource for the population and the natural environment. On the one hand, over the past 
centuries, traditional water use systems (WUSs) such as qanats (a slightly sloping tunnel 
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constructed to accumulate and transfer water from the groundwater to the land surface), wells, 
springs and dams have been used by native residents to manage limited water resources (Motiee et 
al., 2006; Ouessar et al., 2009; Nasiri and Mafakheri, 2015). Qanats and drilled wells are considered 
the most stable and successful water supply systems in hot and arid climates (Boustani, 2008; 
Hussain et al., 2008; Mostafaeipour, 2010). Hence, WUSs are a reliable way to secure access to 
safe and affordable water for several purposes, e.g., drinking, washing and irrigating. On the other 
hand, numerous soil and water conservation measures (SWCMs) such as semi-circular bunds and 
soil bunds have been constructed to decrease the velocity of surface runoff and erosion rate, and to 
consequently increase the infiltration and percolation in the Halilrood Basin of Iran. The existence 
of WUSs and SWCMs in the Wadis has led to complex hydrological impacts that are not well 
understood (Ouessar et al., 2009; Abouabdillah et al., 2014; Hashemi et al., 2015). The WUSs are 
used to not only transfer water from shallow aquifers to the land surface for various purposes, but 
also release water into rivers in some cases, which makes it difficult to include these water usages 
in a hydrological model. Although stream flow may be separately measured for the qanats, springs 
and wells, the amount of released or utilized water from these water supply systems has not been 
recorded so far. Hence, the amount of extracted water that is released to rivers is unknown. In 
addition, in a Wadi catchment, the amount of water which is extracted by these traditional water 
supply systems can significantly affect river runoff (Ouessar et al., 2009). 

In previous studies that have been carried out in Iran, mostly the impacts of WUSs and SWCMs 
on groundwater were evaluated, but they were not included in a catchment model to simulate stream 
flow. For instance, Sadeghi-Tabas et al. (2017) integrated qanats, springs and wells in a 
groundwater model with a genetic multi-algorithm method to define the pumping rates within a 
series of Pareto solutions. Naghibi et al. (2018) used the location of qanats to model the 
groundwater extraction potential of Beheshtabad Basin in the center of Iran. Understanding 
hydrological processes and providing reasonable strategies and plans for policy makers for a better 
management of water resources in these regions are important. In this regard, hydrological models 
can be used to depict complex hydrological conditions and investigate the impacts of WUSs and 
SWCMs on hydrological processes in a Wadi catchment. Due to the specific hydrological and 
climatic conditions and limited data availability, using a model that is capable of simulating 
hydrological processes under these conditions, is critically important. 

Hydrological processes and stream flow of arid and semi-arid regions were simulated by a variety 
of models in different parts of the world, and a large number of studies focused on stream flow 
simulations in dry basins. For instance, Hernandez et al. (2000) applied two hydrological models 
in a small semi-arid watershed in southeastern Arizona, USA; Peugeot et al. (2003) employed the 
r.water.fea model in Sahelian West Niger; Riad et al. (2004) used an artificial neural network (ANN) 
model for the Qurika Wadi Basin in Morocco; McMichael et al. (2006) estimated the monthly 
stream flow of a semi-arid basin in central California using MIKE SHE model; and Lesschen et al. 
(2009) simulated runoff of Carcavo Basin with the LAPSUS (landscape process modelling at multi- 
dimensions and scales) model in Southeast Spain. Also the Soil and Water Assessment Tool 
(SWAT) model was already applied in arid and semi-arid basins with limited data availability in 
different parts of the world (Yebdri et al., 2007; Ning et al., 2015; Cheng et al., 2017; Ignatius and 
Jones, 2017; Hallouz et al., 2018). In addition, the SWAT model has successfully been used in 
different parts of Iran, mostly to simulate stream flow (Rostamian et al., 2008) and groundwater 
(Izady et al., 2015), as well as to estimate sediment (Rostamian et al., 2008) and nitrate transport 
(Jamshidi et al., 2010), and to assess impacts of climate change (Zahabiyoun et al., 2013; Rafiei et 
al., 2015), land use change (Ghaffari et al., 2009; Aghsaei et al., 2020) and engineering projects 
(e.g., dam construction) (Ghobadi et al., 2015) on hydrology. In summary, the SWAT model has 
shown its capability of modeling water fluxes in Iran as well as in arid and semi-arid basins in the 
world. 

Despite considerable progress in hydrological simulations of Wadis during the last two decades 
(Al-Qurashi et al., 2008; McIntyre and Al-Qurashi, 2009; Ben and Abida, 2016), only few efforts 
have been made to include WUSs and SWCMs in hydrological models to simulate stream flow. As 
the SWAT model is an open source model that includes subroutines and parameters related to 
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agricultural management practices as well as water uses, it is often used to consider different WUSs 
and SWCMs in the simulation of stream flow. For instance, Ouessar et al. (2009) applied the SWAT 
model to investigate the effect of WUSs on the water balance components of an arid watershed of 
Wadi Koutine in Southeast Tunisia, based on 38 runoff events between 1973 and 1985. 
Abouabdillah et al. (2014) employed the SWAT model to evaluate the impact of SWCMs (i.e., 
contour ridges) on hydrological components and erosion in Merguellil catchment in the center of 
Tunisia. The model predicted that contour ridges produced an annual reduction in surface runoff, 
an increase in aquifer recharge and a retention of large proportions of entrained sediment. 
Moreover, Khelifa et al. (2017) used the SWAT model to analyze the effect of bench terraces on 
water and sediment yield in an experimental catchment (3.2 km?) and found that the local terraces 
reduced both surface runoff and sediment yield by around 20%. Although the impacts of different 
WUSs and SWCMs on surface runoff are shown in these previous studies, the influences on 
different segments of the hydrograph and on model performance have not been evaluated when 
WUSs and SWCMs are implemented. 

Therefore, the aims of this study are (1) to integrate traditional WUSs and SWCMs in the 
Halilrood Basin of Iran into a hydrological model, (2) to assess the performance of this integration, 
and (3) to quantify and discuss the impacts of this integration on hydrological processes of the 
Halilrood Basin. 


2 Materials and methods 


2.1 Study area 


Halilrood Basin is located in Kerman Province in the center of Iran (Fig. 1) and comprises an area 
of approximately 7224 km?. The Halilrood River is the largest river in terms of discharge in Kerman 
Province and is one of the major water sources for Jazmorian wetland, which is crucial for the 
natural ecosystems in Southeast Iran (Skandari et al., 2016). Due to the mountainous area in the 
north, elevation of the entire basin varies from approximately 1391 to 4359 m a.s.l. About 75% of 
the basin is covered by bare land (ESA, 2010). Limited rainfed agriculture and irrigated farming 
are taking place only in the surroundings of the river and qanat channels. 

Lithosol, Calcaric Regosol and Calcic Yermosols are the dominant soils in the basin, which are 
classified into soil hydrologic group C (FAO/IIASA/ISRIC/ISS-CAS/JRC, 2009) with a slow 
infiltration rate and a high runoff potential. Lithosol and Calcic Yermosols are shallow and 
moderately deep soils, containing a higher sand content and a lower silt and clay content. 
Specifically, the top soil layer consists of a nearly equal percentage of sand (44%) and silt (36%). 
Calcaric Regosol is a deep soil and typically has a clay-loam texture with a high clay content. 

Halilrood Basin is characterized by a desert climate with hot and dry summers according to the 
Köppen-Geiger climate classification. The maximum daily average temperature can reach up to 
40°C at Kenaroyeh station located near the outlet of the basin. The long-term annual average 
temperature is 13°C. The long-term average annual precipitation in the basin is less than 225 mm 
(1993-2009), most of which is received between January and May, whereas precipitation is 
negligible between June and December. The annual potential evapotranspiration (PET) ranges from 
2039 to 2569 mm based on the observation data from the Baft synoptic station (Fig. 1). 

Five hydrometric, nine climatic and one synoptic station are in operation in the upstream of 
Kenaroyeh station, for which daily climate (from 1979 to 2010) and discharge (from 1993 to 2011) 
data are available. The outlet of the basin is located in the upstream of Jiroft dam at Kenaroyeh 
station (shown in Fig. 1). The mean annual discharge from 1993 to 2011 is 7.68 m?/s. In the northern 
part of the basin, Baft and Rabor (Nabi-e-Akram) dams (ratio of reservoir volume to mean annual 
runoff volume is 0.140 and 0.002, respectively) are in operation since 2007 and 2009, respectively 
(IWPCO, 2018). 

In addition to using surface water, extracting groundwater is common to overcome dry periods 
and particularly the drought disasters of the last decade. Hence, the necessary water for domestic, 
industrial and agricultural purposes is supplied from 422 wells, 347 qanats, 2 dams and 184 springs 
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scattered throughout the basin (Fig. 1). 
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Fig. 1 Location of Halilrood Basin in Kerman Province of Iran (a) and water use systems, and climatic and 
hydrometric stations in the Halilrood Basin (b) 


2.2 Hydrological model 


The SWAT model (Arnold et al., 1998) is one of the most widely used catchment models which can 
be applied to simulate hydrological processes in different climatic regions and under various 
conditions. It is a semi-distributed model that splits the basin into several sub-basins, where each 
sub-basin consists of hydrological response units (HRUs) with unique combinations of soil, land 
use and slope (Arnold et al., 2012). The SWAT model has been developed to quantify the impacts 
of water and agricultural management practices on stream flow and other hydrological components. 
The model enables the complex simulation of detailed hydrological process in Wadis through the 
utilization of a wide range of parameters, which however requires experience for a successful 
application of the model. In addition, the model has shown its capability of modeling water fluxes 
in Iran as well as in arid and semi-arid basins elsewhere with the integration of WUSs and SWCMs. 


2.3 Data input and model setup 


2.3.1 Model parameterization 


A QGIS interface for SWAT (QSWAT) model was used to prepare the SWAT model input files (Dile 
et al., 2016). We used the SRTM (shuttle radar topography mission) digital elevation model (Jarvis 
et al., 2008), soil data from the harmonized world soil database (FAO/IIASA/ISRIC/ISS-CAS/JRC, 
2009) and a Globcover 2009 land use map (ESA, 2010) to set up the SWAT model. It should be 
noted that we divided the Halilrood Basin into 285 sub-basins with 6091 HRUs using five slope 
bands (<3%, 3%-5%, 5%-8%, 8%-15% and >15%) based on the FAO classification. This setup 
provided a sufficient spatial resolution for implementing the SWCMs in the SWAT model. 

Precipitation data from nine climatic stations scattered within the Halilrood Basin were used 
(Fig. 1). Additional climate data of wind speed, temperature, solar radiation and relative humidity 
were collected from the Baft synoptic station in the north of the basin (Fig. 1). We determined two 
model setups by running the model without (Default model) and with (WUS-SWCM model) WUSs 
and SWCMs, to analyze the impacts of WUSs and SWCMs in the SWAT model. 

Discharge data from the hydrometric station located at the outlet of the basin (Kenaroyeh station) 
was divided into calibration and validation datasets. We selected the years for calibration and validation 
based on representative climatic conditions in each period, i.e., equal distribution of dry year (total 
precipitation<200 mm), wet year (total precipitation>270 mm) and average year (200-270 mm annual 
precipitation). The calibration and validation periods chosen for the simulation runs were 1995- 
2003 and 2004-2009, respectively, with a two-year spin-up phase prior to 1995. 
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2.3.2 Water use systems (WUSs) 

The WUSs in the Halilrood Basin include qanats, springs, wells and dams. Qanats are regarded as 
a traditional water use system in the study area, and more details can be available in Nasiri and 
Mafakheri (2015) and Mostafaeipour (2010). Qanats, wells and springs are providing water from 
shallow groundwater for domestic and industrial purposes and irrigation in the Halilrood Basin. To 
consider qanats, spring and wells in the SWAT model, we applied consumptive water use (an 
embedded approach in the SWAT model) to remove water from the basin on a monthly time step. 
The average monthly amount of water (m?/d) was therefore removed from the specified source such 
as the shallow aquifer (WUSHAL), the deep aquifer (WUDEEP), the reach (WURCH) or ponds 
(WUPND) within any sub-basins in the basin. All water supply systems (qanats, wells and springs) 
mentioned above were modeled by taking water from the shallow aquifer using the variable 
WUSHAL (Table 1). The actual locations of qanats, wells and springs were used to specify the 
respective sub-basins and HRUs in the SWAT model. It is notable that, due to the structure of 
qanats, the location of the mother well (the last and deepest well) for each qanat was considered to 
find out the origin of the water. Moreover, we summed up the extracted amount of water from the 
shallow aquifer by qanats, wells and springs based on the measured discharge provided by Iran 
Water & Power Resources Development Company (IWPCO) in 2001 and 2006 and included it in 
the SWAT model. Due to changes in groundwater recharge and excessive consumption of 
groundwater, some of the qanats are dry and deactivated nowadays based on the measured 
discharge provided by IWPCO in 2001 and 2006. These were kept out of the calculation (Table 1). 
In addition, the wells were classified into two groups, i.e., shallow and semi-deep wells. Due to a 
drop in groundwater table, 3% of the wells were extended to semi-deep wells based on the measured 
discharge provided by IWPCO in 2001 and these wells were drilled deeper in Quaternary 
sedimentary formations (>80 m). Since those wells extracted water from the deep aquifer, they 
were neglected in the SWAT model. Further information about the implemented qanats, wells and 
springs is given in Table 1. 


Table 1 Characteristics of the water use systems (WUSs) in the study area 


Total Drilling period 
SWAT- Hydrological annual Avage 
Number of f Implemented daily 
WUSs Number : variables components extracted ,. 1966- 1981- 1995- 
deactivated scale discharge 
affected affected water (L/s) 1980 1994 2011 
(x10* m3) 
Qanats 347 85 wusa Groundwater  Sub-basin 16.67 1.52 * * * 
baseflow 
Springs 184 14 WUSHAL Groundwater,  Sub-basin 9.99 1.36 * * * 
baseflow 
Wells 422 93 WUSHAL Groundwater sub-basin 48.56 * 65 59 268 
baseflow 


Note: SWAT, Soil and Water Assessment Tool; WUSHAL, removal water from the shallow aquifer; *, no data available. 


Baft and Rabor dams are two constructed reservoirs in the headwaters of the basin (Fig. 1), which 
supplied water for irrigation and domestic purposes since 2007 and 2009, respectively. Table 2 
shows the reservoir parameters used in the SWAT model. 


2.3.3 Soil and water conservation measures (SWCMs) 


Due to the soil properties and climatic conditions in the Halilrood Basin, semi-circular bunds (so 
called Eyebrows) and soil bunds have been constructed in bare land to protect the soil from erosion 
and to collect surface water. The semi-circular bunds and soil bunds, belonging to SWCMs, can 
alter hydrological processes by reducing surface runoff and by increasing evapotranspiration and 
infiltration. Those might also change the small surface runoff routing schemes. However, since we 
are mainly interested in the impact on the water balance, the influence of SWCMs was estimated 
by changing the curve number (CN) for the affected HRUs, which is in agreement with the studies 
of Arabi et al. (2007), Taye et al. (2013) and Adimassu et al. (2014). 
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Table 2 Hydrological details of the two reservoirs (Baft and Rabor dams) 


Parameter Unit Description Baft Rabor 
SUB-BASIN - Number of the sub-basin in which the reservoir is located 169 198 
TYRES - The operational year 2007 2009 
RES_ESA hm? Reservoir surface area when the reservoir is filled to the emergency spillway 8.0 85.0 
RES_EVOL x10* m? Volume of water needed to fill the reservoir to the emergency spillway 70 4000 
RES_PSA hm? Reservoir surface area when the reservoir is filled to the principal spillway 7.4 32.0 
RES_PVOL x10* m? Volume of water needed to fill the reservoir to the principal spillway 65 3500 
RES_VOL x10* m? Initial reservoir volume 30 2000 


Notes: -, no unit. 


In the SWAT model, SWCMs were implemented through a reduction of the CN values which 
consequently decreased the amount of surface runoff and increased the amount of infiltration. 
Based on previous studies, this surface runoff reduction varies from 28% (Adimassu et al., 2014) 
to 50% (80%) (Taye et al., 2013) for all respective sub-basins and HRUs. In our model, 50% of the 
CN values were reduced in the HRUs with the artificial structures (Table 3). 


Table 3 Characteristics of the soil and water conservation measures (SWCMs) in the study area 


SWCMs Number of sub-basins | SWAT-variables affected Hydrological components and processes affected 
Semi-circular bunds 17 Curve number Surface runoff and infiltration 
Soil bunds 53 Curve number Surface runoff and infiltration 


2.3.4 Calibration 


In this study, we reviewed previously published SWAT model studies in Wadi catchments to depict 
the hydrological components of the model, and carried out a manual sensitivity analysis to select 
the most important hydrological parameters. The properties and variables governing water 
movement into the soil and consequently into or out of the shallow aquifer were the most important 
parameters to sufficiently represent the WUSs in the model. The selected eight parameters and their 
ranges shown in Table 4 were based on previous studies in arid and semi-arid areas (Shrestha et al., 
2016; Qi et al., 2017; Zettam et al., 2017) and the manual sensitivity analysis. 

Latin hypercube sampling (LHS) from the R package FME was used to generate a set of 
variations for the calibration parameters (Soetaert and Petzoldt, 2010; Pfannerstill et al., 2014; Haas 
et al., 2016) for 3000 model runs. The same set of LHS was used for the two setups. For each 
calibration run, the SWAT model input files were rewritten in R software (Pfannerstill et al., 2013; 
R Core Team, 2013). 

Table 4 Selected parameters for calibration in the SWAT model 
Calibration range 


Parameter Description Unit = > Type 
Minimum Maximum 


CN? Initial soil conservation service runoff curve z -30 “15 Add 
number for moisture condition II 


mm H,O/mm 


SOL_AWC Available water capacity of soil layer soil —0.5 0.5 Add 
ESCO Soil evaporation compensation factor - 0.90 0.96 Replace 
GW_DELAY Ground water delay time d 4 10 Replace 
RCHRG_DP Deep aquifer percolation fraction - 0.5 0.9 Replace 
ALPHA_BF Base flow alpha factor Per day 0.08 0.20 Replace 

SOL_K Saturated hydraulic conductivity mm/h 30 40 Add 
EVRCH Reach evaporation adjustment factor - 0.5 0.8 Replace 


Notes: -, no unit. Add means that a given value is added to the existing parameter value; Replace means that the existing parameter value 
is to be replaced by a given value. 


2.3.5 Model evaluation 
The Nash-Sutcliffe efficiency (NSE; Nash and Sutcliffe, 1970), percent bias (PBIAS; Gupta et al., 
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1999), root mean squared error (RMSE) and the ratio of standard deviation (RSR; Moriasi et al., 
2007) are several quantitative criteria that are frequently used to evaluate the performance of 
hydrological models. In addition, the Kling—Gupta efficiency (KGE) is a statistical performance 
metric that considers bias, correlation and variability separately (Gupta et al., 2009). Kling et al. 
(2012) published a modified version of the KGE, in which bias and variability ratios were not cross- 
correlated. 

Since each of the aforementioned criteria has a specific hydrological focus (Guse et al., 2019), 
SWAT model parameters were calibrated using a multi-metric approach which has been proven 
efficient to balance model performance (Pfannerstill et al., 2014; Haas et al., 2016; Tigabu et al., 
2019). Accordingly, 3000 model runs were carried out with the SWAT model. To assess model 
performance, we used NSE, PBIAS, RSR and the modified KGE as the performance measures on 
a daily basis. To identify the best calibration runs for both model setups and to enable a comparison, 
we defined thresholds for NSE, PBIAS and RSR so that at least 5% of the total model runs were 
remaining. To this end, a hierarchical selection of model runs was conducted. Firstly, the model 
runs with NSE values greater than or equal to 0.68 were selected. Secondly, the model runs with 
PBIAS between —25% and 25% were selected. Thirdly, the model runs with RSR greater than or 
equal to 0.56 were identified. After the application of these thresholds on the 3000 model runs, we 
sorted the selected model runs according to the KGE. Finally, we selected 5% of the total model 
runs (150 model runs) as the best calibration runs based on the KGE for both model setups. 

In this study, we used the segmentation of the flow duration curve (FDC; Yilmaz et al., 2008) 
and performance criteria for every segment (Pfannerstill et al., 2014) to distinguish the impact of 
WUSs and SWCMs on different parts of the hydrograph. The segments were split at different 
exceedance probabilities of the FDC: 0-5%, 5%-20%, 20%-—70%, 70%-95% and 95%-—100%, 
which were associated to very high, high, middle, low and very low flows, respectively. Equal 
ranges were considered for the very low and very high flows in the FDC as described in Pfannerstill 
et al. (2014). The RSR was applied on each of the five segments, which enabled a tailored 
evaluation of the model performance for the 150 best calibration runs (Haas et al., 2016). Since the 
WUSs (qanats, wells and springs) extract water from the shallow aquifer and SWCMs increase the 
infiltration rate, impacts will likely be highest on the base flow. Therefore, RSR in combination 
with different parts of the FDC is a suitable measure to evaluate the effects of WUSs and SWCMs 
on the model performance. 


2.3.6 Water balance components 


To assess via which pathways stream flow is affected by implementing WUSs and SWCMs in the 
model, we compared the hydrological components water yield (WYLD) and actual 
evapotranspiration (ET) of the 150 best calibration runs. WYLD represents the total amount of 
water leaving the sub-basin and entering the main channel, which can be evaluated by the model 
as follows (Neitsch et al., 2011): 

WYLD=SURFQ+LATQ+GWQ-TLOSS, (1) 
where WYLD is the water yield (mm); SURFQ is the surface runoff (mm); LATQ is the lateral flow 
contribution to stream (mm); GWQ is the groundwater contribution to stream flow (mm); and TLOSS 
is the transmission losses (mm), i.e., water loss via transmission through the bed of the channels. 
Since the parameter combinations of the best calibration runs may not be the same, the parameter 
combinations of the best calibration runs were used for both model setups separately and the average 
changes for the aforementioned hydrological components were compared for all the selected runs. 


3 Results and discussion 


3.1 Calibration 


Figure 2 shows the results of the performance metrics selection of both model setups for the 
calibration period of 1995-2003. In Figure 2, the 3000 model runs are shown as a relationship 
between each performance metric and KGE as the main performance criteria. For each performance 
metric, the acceptable calibration runs are identified. For the Default model, the KGE values of the 
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Fig.2 Number (n) of acceptable calibration runs for each performance metric (black points) for the Default model 
in the left column and the WUS-SWCM model in the right column. Default model, a model setup without water 
use systems and water conservation measures; WUS-SWCM model, a model setup with water use systems and 
water conservation measures. The gray points represent the range of the KGE (Kling—Gupta efficiency) for the 
complete dataset of the 3000 model runs. The last row of both columns shows the selection of acceptable calibration 
runs after the application of all thresholds for the different performance metrics. NSE, Nash-Sutcliffe efficiency; 
PBIAS, percent bias; RSR, ratio of standard deviation. 


3000 model runs range from —0.1 to 0.8 as compared to a narrower range of the KGE values for 
the WUS-SWCM model, in which the KGE values are mostly greater than 0.3 and 97% of the 
model runs have KGE values higher than 0.5. The number of acceptable calibration runs (n) for 
each metric is defined and highlighted in black for both model setups. For the Default model, the 
number of acceptable calibration runs with NSE>0.68 is 243. The application of the PBIAS 
threshold reduces the number of total model runs from 3000 to 2358 acceptable calibration runs. 
Moreover, the number of acceptable calibration runs is even smaller (182) when the defined 
threshold of RSR is set. However, for the WUS-SWCM model, a higher number of simulation runs 
remains for each performance metric when the threshold is applied. The number of acceptable 
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calibration runs with NSE>0.68 is 2486. The PBIAS threshold leads to a high number of acceptable 
calibration runs (2437). When applying the RSR threshold, only 737 runs are rejected. 

Since the thresholds were applied based on performance metrics of the multi-metric framework 
(NSE, PBIAS and RSR), 177 model runs are remaining for the Default model, and this number is 
higher for the WUS-SWCM model, with 1942 model runs. 

The acceptable calibration runs of the Default model and WUS-SWCM model setups are sorted 
in increasing order of the KGE values for both model setups and finally 150 best calibration runs 
are selected and plotted in Figure 3. The curves of the KGE clearly show the characteristics of the 
two model setups. In both model setups, a good performance is achieved for the first 20 runs, in 
which the KGE value in the calibration period is approximately 0.9. However, the KGE value of 
the 150" Default model run drops below 0.7. In contrast, the minimum KGE value is higher than 
0.8, even for the 150" WUS-SWCM model run (Fig. 3). 


— Calibration (WUS-SWCM model) -—Validation (WUS-SWCM model) 
Calibration (Default model) Validation (Default model) 


0.9 


KGE 


0.4 
0 10 20 30 40 50 60 70 80 90 100 110 120 130140 150 
Number of the best calibration runs 


Fig. 3. Comparison of the 150 best calibration runs for the KGE values of the Default model and the WUS-SWCM 
model setups in the calibration (solid lines) and validation (dashed lines) periods 


Although the KGE values are lower in the validation period when compared to the calibration 
period, a similar pattern is found for the validation and calibration periods (Fig. 3). While the first 
eight runs of the Default model setup reach higher KGE values, the KGE values of the Default 
model setup declines more rapidly as compared to the WUS-SWCM model setup. In general, Figure 
3 shows that better model parameterizations are found for the model with WUSs and SWCMs than 
for the one without. 

The results show that the overall KGE significantly increases when WUSs and SWCMs are 
included. To show the effect of WUSs and SWCMs on different segments of the hydrograph and 
possible improvement, we plotted the ranges of the flow duration curves (FDCs) of the selected 
150 best calibration runs for both model setups together with the observed FDC (Fig. 4). We also 
calculated and summarized the average RSR value over the 150 best calibration runs to quantify 
the difference of both model setups for each FDC segment. The results are shown in Table 5. For 
the calibration period, although a high goodness of fit is achieved for very high flow for both model 
setups (Fig. 4b), the Default model performs better with an average RSR value of 0.22. High flow 
is underestimated by the Default model. This is improved when the WUSs and SWCMs are 
included in the SWAT model (Fig. 4c), which is confirmed by the comparison of the average RSR 
values for the high flow segment in Table 5, where RSR values decrease from 0.42 to 0.22. Flow 
in the middle segment is overestimated when the WUSs and SWCMs are implemented in the SWAT 
model (Fig. 4d), which results in higher RSR values. The average RSR value decreases from 3.10 
to 2.60 (Table 5) for the low flow segment. 

In contrast, for the validation period, when the WUSs and SWCMs are implemented, the model 
is improved for very high and middle flows (Figs. 4h and j), with RSR values reducing from 0.32 
to 0.28 and from 0.87 to 0.81, respectively (Table 5). The comparison of the average RSR value for 
the high flow segment reveals a worse performance for the WUS-SWCM model (Table 5), with an 
overestimation in the high flow segment (Fig. 4i). 
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Fig. 4 Flow duration curves (FDCs) of the selected 150 best calibration runs for the WUS-SWCM model (light 
green) and Default model (pink) in the calibration (1995-2003; a-f) and validation (2004-2009; g-l) periods. The 
common FDCs between the two model setups are shown in dark green. Different segments of the hydrograph are 
separated by dotted vertical lines. 


In both model setups, the model is not capable of estimating low and very low flows close to 
zero (Figs. 4e, f, k and 1). This means that low and very low flows are simulated as zero, while 
these values are only close to zero in the observed data. This leads to the same RSR value for very 
low flow segment in both calibration and validation periods. Some qanats have a pool at the end of 
the channel that allows for a slower release of the water. As this has not been included in the model 
implementation of qanats, which might explain the incoherencies in low and very low flows. 

The parameters and their values which lead to the best model performance for both model setups 
are represented in Table 6. While the values of some parameters such as the soil evaporation 
compensation factor (ESCO) and deep aquifer percolation fraction (RCHRG_DP) are almost the 
same in both model setups, other parameters differ considerably. For instance, CN value and reach 
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evaporation adjustment factor (EVRCH) are higher in the WUS-SWCM model, i.e., more surface 
runoff and less infiltration as well as more evaporation from the reach are possible in the WUS- 
SWCM model. In those areas where SWCMs are present, the CN value is further reduced in the 
WUS-SWCM model setup prior to calibration. This explains the slightly lower overall reduction 
in the WUS-SWCM model setup compared to the Default model during calibration, as shown by 
the parameter set. Moreover, groundwater flows are faster in the WUS-SWCM model with a higher 
ground water delay time (GW_DELAY) value. 


Table 5 Summary of the application of the ratio of standard deviation (RSR) for each flow duration curve (FDC) 
segment for the average of the 150 best calibration runs of each model setup in calibration (1995-2003) and 
validation (2004—2009) periods 


Calibration period (1995-2003) Validation period (2004—2009) 
Performance metrics Performance metrics 
RSR RSR 
KGE 
Very high High Middle Low Very low Very high High Middle Low Very low 
Default model 0.77 0.22 0.42 0.38 3.10 3.71 0.58 0.32 0.53 0.87 3.44 0.99 

WUS-SWCM model 0.82 0.29 0.22 0.70 2.60 Skil 0.61 0.28 1.04 0.81 3.44 0.99 

Relative changes te — + — sr ve ar te — qP z $ 


Note: Default model, a model setup without WUSs and SWCMs; WUS-SWCM model, a model setup with WUSs and SWCMs. + means 
a better performance in the WUS-SWCM model; — means a worse performance in the WUS-SWCM model; * means no change in the 
model performance. 


Model setup 
KGE 


Table6 Parameter sets that lead to the best model performance for each model setup 


Number of Parameters 
model runs CN2 SOL_AWC ESCO GW_DELAY RCHRG_DP ALPHA_BF SOL_K EVRCH 


(add) (add) (replace) (replace) (replace) (replace) (add) (replace) 
Default model 292 -24.610 0.005 0.956 8.345 0.527 0.143 32.359 0.655 


WUS-SWCM model 2868 -20.917 —0.008 0.954 4.791 0.508 0.191 36.070 0.798 


Model setup 


To understand how the WUSs and SWCMs affect the hydrograph, we compared the hydrographs 
of the two model setups using the parameter sets of the best WUS-SWCM model (Fig. 5). Stream 
flow and particularly stream flow peaks are overestimated by the Default model (Fig. 5a). We 
compared the observed and simulated stream flow for a shorter period (January—June of 2001) from 
both model setups to assess differences between the two model setups in more detail (Fig. 6). The 
Default model overestimates the peak flows considerably, while implementing WUSs and SWCMs 
results in decreased peak flows and a more reasonable performance in mid-January and at the end 
of February. Low flow is considerably underestimated by the Default model and a higher goodness 
of fit is obtained when the WUSs and SWCMs are implemented, e.g., in February and March. Also, 
the falling limb of the hydrograph is better represented by the WUS-SWCM model. This might be 
due to the higher infiltration rate and a slower release of water to the reach in the WUS-SWCM 
model. From the beginning of the dry period (May), both model setups perform similarly. 


3.2 Hydrological components 


The comparison of hydrological components of the model setups shows that aggregated basin 
value changes are small (changes in ET and WYLD are lower than 3%), because they probably 
compensate for each other on the large scale. Therefore, the change of hydrological components is 
evaluated at the sub-basin scale by comparing the WUS-SWCM model with the Default model 
(Fig. 7a). WYLD is decreasing in all sub-basins with WUSs and SWCMs, varying from 0 to -18% 
(Fig. 7b). The change in runoff components is split in surface runoff and groundwater flow (Figs. 
7c and d). It should be noted that evapotranspiration and lateral flow changes vary less than 1% 
and therefore are not shown in this paper. The surface runoff ranging from —30% to —99% and its 
contribution to the stream decrease in sub-basins with WUSs and SWCMs. Moreover, a higher 
contribution of groundwater to the streams in most sub-basins with SWCMs and on the contrary a 
lower groundwater contribution in sub-basins with WUSs (qanats, wells and springs) are shown, 
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Fig. 5 Comparison of observed and simulated stream flow from the Default model (a) and the WUS-SWCM 
model (b) setups in the calibration (1995-2003) and validation (2004—2009) periods 
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Fig.6 Detailed comparison of observed and simulated stream flow from the Default model and the WUS-SWCM 
model setups for the period of January-June in 2001 


but in sub-basins with both WUSs and SWCMs, groundwater contribution increases (maximum 
46%) or does not show any change. This indicates that SWCMs counterbalance the impact of WUSs 
on groundwater flows. Although WUSs are taking water from the shallow aquifer and change 
groundwater flows, the implementation of SWCMs in the basin compensates for the extracted water 
from the shallow aquifer by decreasing the surface runoff and increasing the infiltration rate and 
ground water recharge. This finding is in agreement with the studies of Abouabdillah et al. (2014) 
and Khelifa et al. (2017). 


4 Conclusions 


In this study, WUSs and SWCMs were successfully implemented into a hydrological model in the 
Halilrood Basin. Our results clearly show that there are differences in both model setups (Default 
model and WUS-SWCM model), which can be related to the implementation of the different 
measures. Model performance improved when WUSs and SWCMs are included in the SWAT 
model. However, this improvement is not similarly observed in all segments of the hydrograph. 
The model is capable of simulating the hydrological processes more realistically when more details 
of water systems are considered. The comparison of the hydrological components illustrates that 
the contribution of surface runoff to the stream is decreased pronouncedly when WUSs and 
SWCMs are implemented. The presentation of SWCMs leads to a higher contribution of 
groundwater to the stream flow by increasing the infiltration rate and groundwater recharge and 
decreasing the surface runoff. Therefore, in the sub-basins with both WUSs and SWCMs, the 
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Fig. 7 Spatial distribution of sub-basins with or without WUSs and SWCMs (a) and changes of hydrological 
components when WUSs and SWCMs are implemented (b-d). (b), water yield; (c), groundwater flow; (d), surface 
runoff. Decreasing change is shown as negative value (red) while increasing change is shown as positive value 
(blue). 


implementation of SWCMs compensates for the impact of the WUSs. Furthermore, the amount of 
water leaving the basin is decreased due to the implementation of WUSs and SWCMs in the model, 
leading to a higher water use in the basin. Considering WUSs and SWCMs in the model is 
recommended to accurately simulate the hydrological processes, particularly in Wadis, where water 
availability is limited, these measures have a pronounced effect. Moreover, with regard to scenario 
assessments in the arid and semi-arid regions, hydrological models that include these measures are 
needed to yield reliable results and for being able to evaluate the actual impact of climate change. 
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